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ABSTRACT 

We report on VERITAS very-high-energy (VHE; E>100 GeV) observations 
of six blazars selected from the Fermi Large Area Telescope First Source Catalog 
(1FGL). The gamma-ray emission from 1FGL sources was extrapolated up to the 
VHE band, taking gamma-ray absorption by the extragalactic background light 
into account. This allowed the selection of six bright, hard-spectrum blazars 
that were good candidate TeV emitters. Spectroscopic redshift measurements 
were attempted with the Keck Telescope for the targets without Sloan Digital 
Sky Survey (SDSS) spectroscopic data. No VHE emission is detected during 
the observations of the six sources described here. Corresponding TeV upper 
limits are presented, along with contemporaneous Fermi observations and non- 
concurrent Swift UVOT and XRT data. The blazar broadband spectral energy 
distributions (SEDs) are assembled and modeled with a single-zone synchrotron 
self-Compton model. The SED built for each of the six blazars show a synchrotron 
peak bordering between the intermediate- and high-spectrum-peak classifications, 
with four of the six resulting in particle-dominated emission regions. 

Subject headings: active galactic nuclei: general — gamma rays: individual(RGB 
J0136+361, RGB J0316+090, RGB J0909+231, RGB J1058+564, RGB J1243+364, 
RX J1436. 9+5639) 
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1. Introduction 

Blazars are active galactic nuclei (AGN) with a relativistic jet pointed close to the Earth 
line of sight. AGN are thought to be powered by accretion onto a supermassive black hole 
(~ 1O 9 M 0 ) at the center of the host galaxy and are characterized by a double-peaked spectral 
energy distribution (SED) in the i'Fj, representation. 

The lower energy peak of the broadband SED is attributed to synchrotron emission of 
highly relativistic electrons and positrons in the presence of a tangled magnetic field. In lep- 
tonic models, the higher energy peak is produced via the inverse-Compton (IC) up-scattering 
by the relativistic leptons of the synchrotron photons (synchrotron self-Compton, SSC) or 
a photon field external to the jet (external Compton, EC). More details regarding leptonic 
non-thermal emission of blazars can be found in Dernier et al. (1992); Maraschi et al. (1992); 

Sikora et al. (1994) and the references therein. Alternative models associate the higher- 
energy peak to interaction of relativistic protons with an ambient photon held (Aharonian et al. 

2002; Bednarek 1993; Dar et al. 1997; Mannheim 2000; Miicke & Protheroe 2000; Pohl & Schlickeiser 
2000), or a hybrid population comprised of both leptons and hadrons (Bottcher 2007). 

The blazar population is divided into two subclasses: hat spectrum radio quasars (FS- 
RQs) and BL Lacs. FSRQs are, to hrst order, more distant, more luminous, and have 
stronger emission lines than BL Lacs. Historically, BL Lacs have been sub-classihed based 
on their radio and X-ray hux ratios as low-, intermediate- or high-frequency-peaked BL Lac 
objects (LBL, IBL and HBLs, respectively) as described in Padovani & Giomrni (1995) and 
Bondi et al. (2001). More recently, a classification based on the location of the synchrotron 
peak in frequency space {y sy nch) has been proposed by Abdo et al. (2010a), with low- 
spectrum-peaked (LSP) BL Lac objects having v sync h below 10 14 Hz, intermediate-spectrum- 
peaked (ISP) BL Lac objects peaking between 10 14 and 10 15 Hz and high-spectrum-peaked 
(HSP) BL Lac objects showing a peak above 10 15 Hz. HSP BL Lacs are the most common 
extragalactic object to be detected at very-high-energies (VHE; E > 100 GeV), compris- 
ing 33 of 41 VHE BL Lac objects detected as of June 2012. There have also been 4 ISPs 
(Acciari et al. 2009a, b, 2010; Ong 2010a) and 4 LSPs (Mazin et ah 2009; Mariotti 2011; 

Hoffman 2010; Albert et al. 2007) detected since the advent of VHE gamma-ray astronomy 
in the late 1960s. In addition, three FSRQs, three radio galaxies, and two starburst galaxies 
complete the catalog of associated extragalactic objects (Horan & Wakely 2008) 1 . 

The potential scientihc impact of increasing the catalog of VHE emitting extragalactic 
objects is substantial. A significant fraction of the power released from these objects is within 


1 http : //teveat . uchicago . edu/ 
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the VHE band. This makes the measurement of VHE blazar spectra an important component 
of the overall understanding of these objects. With a better sample of well-measured VHE 
blazar spectra available for study, a population-based investigation of gamma-ray production 
in these objects through broadband SED modeling will be possible, providing means to 
answer the long-standing question of whether VHE gamma-ray emission results from leptonic 
or hadronic processes in AGN jets. We can also apply the model-inferred properties of these 
new discoveries to understanding how the gamma-ray production may differ among blazar 
subclasses, exploring the apparent blazar sequence and evolution of the AGN (Fossati et ah 
1998; Ghisellini & Tavecchio 2008; Meyer et ah 2011). 

VHE blazars can also be used to constrain the optical to near-infrared extragalactic 
background light (EBL) density and evolution, as well as the nature of cosmic rays. The 
EBL encodes the integrated history of structure formation and the evolution of stars and 
galaxies in the Universe. Understanding these characteristics requires detailed theoretical 
modeling of all the processes that contribute, e.g. structure formation and stellar evo- 
lution (Dominguez et ah 2011; Gilmore et ah 2009; Primack et ah 2005; Finke et ah 2010; 
Franceschini et ah 2008; Hauser & Dwek 2001; Stecker et ah 2006). According to some cos- 
mic ray models, e.g. Essey & Warren (2012), interactions of cosmic rays along the blazar 
line-of-sight can produce relatively hard VHE gamma-ray spectra as compared with the 
high-energy gamma-ray spectra, depending on the distance to the blazar. 

The current catalog of extragalactic VHE objects largely contains relatively nearby 
AGN; only three have a confirmed redshift above z = 0.3. The highest redshift blazar de- 
tected at VHE thus far is the FSRQ 3C279, at a redshift of 0.536 (Albert et ah 2008b). The 
proximity of these VHE blazars is partly a result of gamma-ray absorption by the EBL. VHE 
gamma rays that propagate through the intergalactic medium are absorbed by low energy 
EBL photons via pair production, 7 + 7 — * e + + e~ (Nikishov 1962; Gould & Shreder 1967; 
Stecker et ah 1992). The absorption process deforms the VHE gamma-ray spectra emitted 
by extragalactic objects in a redshift-dependent manner and can be translated to upper 
limits on the local density of the far-infrared EBL (Stecker et ah 1993; Dwek & Krennrich 
2005; Aharonian et ah 2006; Primack et ah 2011; Orr et ah 2011). 

The search for new VHE blazars is complicated by the fact that many of these objects 
do not yet have known redshifts. BL Lac objects, by definition, display very weak or no 
optical emission or absorption lines used for spectroscopic redshift measurements. Due to 
the interaction of VHE photons with the EBL, non-detection of blazars with no redshift 
information can be either attributed to the object being too distant or having an intrinsically 
low luminosity in the VHE band. This makes VHE discovery observations of blazars with 
no measured redshift a risky venture, although it has proven successful in the past, as 
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in the discoveries of VHE emission from 3C66A, PKS 1424+240, and RX J0648. 7+1516 
(Acciari et al. 2009a, 2010; Aliu et al. 2011). 

The small field of view of imaging atmospheric Cherenkov telescopes (IACTs, < 5.0°) 
makes new source discovery from a large-scale sky survey difficult and therefore the hunt for 
VHE emitting objects has historically involved targeted observations of source candidates 
selected from surveys at lower frequencies, such as the selection of hard X-ray candidates 
presented in Costamante & Ghiscllini (2002). In this way, VHE blazar candidate selection 
has relied on experiments such as EGRET onboard the Compton Gamma Ray Observatory 
(Thompson et al. 1993), Swift (Gehrels et al. 2004) and ROSAT (Turriziani et al. 2007). 
The launch of Fermi in June 2008 has enhanced VHE blazar discovery programs, leading 
to new blazar discoveries including RBS0413, RX J0648. 7+1516, 1ES 0033+595 and 1RXS 
J101015. 9-311909 (Aliu et al. 2012, 2011; Mariotti 2011; Abramowski et al. 2012). 

The Fermi Large Area Telescope (LAT) observes the entire sky in the energy range 
from 20 MeV to > 300 GeV every three hours and has better sensitivity than its predecessor 
EGRET. Within the first 11 months of operation, the 1FGL catalog reported the detection 
of 1451 sources at a significance greater than 5a (Abdo et al. 2010b). A majority of these 
sources are, or are expected to be, associated with AGN (Abdo et al. 2010c). It is known 
that a large majority of these blazars has not yet been detected by VHE instruments, as can 
be seen by the current TeV catalog, which only contains 51 associated extragalactic objects. 
However, the proximity of the Fermi - LAT energy band to that of the IACTs makes the 
1FGL catalog a good place to search for candidate VHE blazars. We report on VERITAS 
observations above 100 GeV of six candidate blazars. For the first time, a multiwavelength 
description of their SED including radio, UV, X-ray and gamma-ray frequencies is assembled 
and modeled using a SSC model. 


2. Target Selection 

The energy coverage of IACTs overlaps with that of Fermi - LAT above 100 GeV and 
extends to tens of TeV. Sources that are most likely to be detected by both Fermi - LAT 
and IACTs have high fluxes and hard photon indices in the Fermi - LAT energy band. A 
selection process was established to identify the best VHE candidates within the first Fermi- 
LAT catalog. Sources at low Galactic latitude (|b| < 10°) were excluded, with the intent of 
removing the majority of Galactic sources from the selection. Additionally, sources with low 
integrated flux above 100 MeV (F 100 MeV < 2 x 10” 9 ph cnrV 1 ), with soft photon indices 
(T > 2.0), or with a low number of associated photons ( N pred < 20), were excluded from the 
selection. Just over 200 sources passed this initial set of cuts. 
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The remaining sources were re-analyzed using the same data set as that used for the 
1FGL catalog to search for the presence of curvature in their spectra. A log-parabola pa- 
rameterization for the spectra was chosen as the alternative to the power-law model (null 
hypothesis). In addition, the data were analyzed only using data above 1 GeV to confirm 
that the results from the power-law fit found for the entire Fermi - LAT energy band agreed 
well with the power-law fit above 1 GeV. Sources showing significant spectral curvature or 
softening of the spectrum above 1 GeV with an improvement in likelihood value correspond- 
ing to 3(7 were excluded from the final selection as these sources are not expected to exhibit 
bright TeV flux levels. For sources with spectra best fit by a power law that also matched 
the power-law fit above 1 GeV, the Ferrm-LAT spectrum was extrapolated up to the VHE 
regime (150 GeV - 1 TeV). When a reliable redshift measurement was available, the EBL 
model of Franceschini et al. (2008) was used to estimate the extrapolated VHE gamma-ray 
flux. Otherwise an assumed redshift of z = 0.2 was used, a conservative value considering 
that most known TeV sources have redshifts less than z = 0.2. Finally, the sources were 
ranked based on their extrapolated flux in the VHE regime. Six VHE candidates, all BL Lac 
objects, were selected for observations with the VERITAS telescope between September 2009 
and June 2010 for 10 hours each, corresponding to the VERITAS 3% Crab flux sensitivity 
exposure timescale. 

Table 1 shows the VHE extrapolated integral flux from the 1FGL catalog power-law 
fits used for selection of the six candidate VHE emitting blazars. These relatively high 
integral flux values above 150 GeV, shown in percentages of the Crab Nebula flux 2 , are 
shown in comparison to the updated power-law fits from 29 months of Fermi-LAT data for 
the steady sources, and data from the VERITAS coincidental window for variable sources. 
The analysis of the extended Fermi-L AT dataset (2008 August 4 to 2011 January 4) is 
detailed in section 4.2. These updated extrapolations show much lower expected integral 
flux values for each of the six candidates, reflecting the fact that as more LAT data were 
collected, better high-energy statistics showed the candidates to be softer and/or dimmer 
than found with the elevent month data set used in the 1FGL. Additionally, we compare 
these extrapolated values to the upper limits derived from the VERITAS observations, where 
the analysis leading to these upper limits is detailed in Section 4.1. 


2 F1ux calculated according to the curved power law presented in Albert et al. (2008a) 


3. The Targets 


RGB J0136+391 (1FGL J0136. 5+3905) was discovered in the third Bologna sky sur- 
vey of 408 MHz radio objects (Ficarra et al. 1985). It was later detected as a X-ray bright 
active galaxy in the Northern ROSAT all-sky survey (Brinkmann et ah 1997), and identified 
spectroscopically as a BL Lac object, with an IBL sub-classification from Laurent-Muehleisen et ah 
(1998) and Laurent-Muehleisen et al. (1999), respectively. Based on the optical and X-ray 
spectral properties, this blazar was proposed as a VHE candidate blazar by Costamante & Ghisellini 
(2002), under the assumption that the redshift was less than 0.2. Bright gamma-ray emis- 
sion above 1 GeV was detected from this source by Fermi after three months of operation. 

Only an upper limit below 1 GeV was reported (Abdo et al. 2009). This blazar also showed 
constant emission in the first eleven months of LAT operation (Abdo et al. 2010b). This 
blazar remains without a known redshift, with no previous spectroscopic redshift measure- 
ments found in the available literature. We attempted a spectroscopic redshift measurement 
on 2009 September 17 (MJD 55091) using the Keck Low Resolution Imaging Spectrometer 
(LRIS) instrument, but measured only a featureless power-law spectrum that is characteris- 
tic of BL Lac objects (see Figure 1 upper-left panel). Detailed spectroscopic analysis of this 
data can be found in Kaplan et al. (2012). 

RGB J0316+090 (1FGL J0316. 1+0904) was first detected by the NRAO Green Bank 
91 meter radio telescope (Becker et al. 1991). It was later optically identified as a BL Lac 
object (Fischer et al. 1998) and sub-classified as an IBL by Laurent-Muehleisen et al. (1999). 

No spectroscopic redshift measurements had been made prior to this work. We attempted a 
spectroscopic redshift measurement on 2011 March 5 using the Keck Echelette Spectrograph 
and Imager (ESI) instrument, resulting in a featureless power-law spectrum except for two 
unidentified absorption features (see Figure 1 upper-right panel). The detailed spectroscopic 
analysis of the Keck ESI data can be found in Kaplan et al. (2012). 

RGB J0909+231 (1FGL J0909. 2+2310) was first detected by the NRAO Green Bank 
radio telescope (Becker et al. 1991) and was later classified as a radio-loud active galaxy 
(Brinkmann et al. 1997). The BL Lac optical counterpart was identified nearly a decade 
later (Mickaelian et al. 2006). The redshift reported by the NASA Extragalactic Database 
(NED) 3 could not be found within the corresponding reference (z = 0.231, Brinkmann et al. 
(2000)) and was therefore taken as unknown. Inspection of publicly available Sloan Digital 
Sky Survey (SDSS) data revealed two Mg II absorption lines in the otherwise featureless 
optical spectrum. Assuming that these lines could be intrinsic, or due to the absorption by 
an intervening cloud, a lower limit on the redshift of z > 0.43 is derived (see Figure 1 rniddle- 


3 http : //nedwww . ipac . caltech.edu/ 
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left panel). This is the only source out of the six selected which has a neighboring Fermi- 
LAT detected blazar within the VERITAS 3° field of view. 2FGL J0910. 9+2246, associated 
with TXS 0907+230, is located 0.61° away from RGB J0909+231. The high redshift blazar 
TXS 0907+230 (z = 2.66 according to Healey et al. (2008)) shows no signal in the 10- 
100 GeV band, with an upper limit of 9.6xl0“ n photons cm~ 2 s _1 reported in the 2FGL 
catalog (Nolan et al. 2012). An exclusion region of radius 0.3° centered on the blazar was 
nevertheless used in the VERITAS analysis to avoid possible contamination. 

RGBJ1058+564 (1FGL J1058. 6+5628) was first detected in the 6C radio survey 
(Hales et al. 1990). It was identified as a BL Lac object during the association of the ROSAT 
all-sky survey with the Hamburg Quasar Survey (Nass et al. 1996). Bondi et al. (2001) clas- 
sified the object as an IBL based on the optical and X-ray characteristics. This blazar was 
detected within the first three months of Fermi - LAT operation (Abdo et al. 2009) with a 
broadband SED shown in Abdo et al. (2010a) which does not include any VHE informa- 
tion. The blazar also shows a high level of flux variability (probability of variability: 79%) 
in the high-energy gamma-ray band, as shown in Abdo et al. (2010b) and is the only one 
out of the six targets in this paper that has a redshift measurement. The absorption lines 
corresponding to the redshift 0.143 can be seen in the SDSS spectrum shown in Figure 1 
(middle-right panel). 

RGB J1243+364 (1FGL J1243. 1+3627) was first reported in the B2 catalog of ra- 
dio sources (Colla et al. 1973). This target was also determined to be a radio-loud ac- 
tive galaxy by ROSAT (Brinkmann et al. 1997), and specifically classified as a BL Lac in 
Appenzeller et al. (1998). NED cited the SDSS data for a spectroscopic redshift of z = 1.065. 
Inspection of this publicly available SDSS data revealed no lines suggesting this redshift but 
instead revealed Mg II absorption lines that translate to a lower limit of z > 0.485 (see Fig- 
ure 1 lower-left panel). A recent photometric redshift of z= 0.51 q;]J from Meisner & Romani 
(2010) is in agreement with the SDSS lower limit. 

RX J1436. 9+5639 (1FGL J1437. 0+5640) was detected in the ROSAT all-sky sur- 
vey and identified as a BL Lac object by Nass et al. (1996). This HBL, as classified by 
Nieppola et al. (2006) based on the frequency of the synchrotron peak, remains without a 
redshift, although the redshift of the galaxy cluster within the same region of the sky is 
known to be z= 0.15 (Bauer et al. 2000). Inspection of the publicly available SDSS data 
shows a featureless spectrum, shown in Figure 1. 
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4. Multiwavelength Observations and Analysis 
4.1. VERITAS 

The VERITAS observatory is an array of four 12-meter diameter IACTs, located in 
southern Arizona. VERITAS is sensitive to photons between 100 GeV and several tens of 
TeV with an energy resolution of better than 20%. The instrument has a 5a point source 
sensitivity of 1% of the Crab Nebula flux in less than 30 hours with an angular resolution of 
less than 0.1° for a Crab- like source with a spectral index of 2.5. See Weekes et al. (2002) 
and Holder et al. (2006) for a detailed overview of the instrument. 

The VERITAS observations of the six VHE candidate blazars were completed between 
September 2009 and June 2010 (MJD 55122-55383). These observations were taken in wobble 
mode, with an offset of 0.5° from the source position in each of four cardinal directions 
to allow simultaneous background measurement, as explained in Fomin et al. (1994) and 
Berge et al. (2007). The radio location of the counterparts as specified by NED were used 
for source position. 

Air shower events initiated by gamma and cosmic-rays are reconstructed following the 
procedure outlined in Acciari et al. (2008). The recorded shower images are parameterized by 
their principal moments, giving an efficient method for suppression of the far more abundant 
cosmic-ray background. Any events with a total charge less than 50 photoelectrons are 
removed from the analysis. A set of cuts is then applied to the parameters in order to 
reject background events (see details of this method in Krawczynski et al. (2006)). These 
parameters reject cosmic-ray like events having the mean scaled width and length of the event 
camera image smaller than 1.1 or 1.4, respectively. Additionally, the reconstructed altitude 
of the maximum Cherenkov emission is required to be higher than 8 km above ground level. 
Gamma-ray like events are extracted from a signal region with a radius of 0.14°, centered at 
the coordinates of the candidate source. 

The quality-selected livetime collected for each target ranges from 4.4 to 14.2 hours 
and results in no detections, with significances ranging from -1.1 to 0.9a, calculated with 
Equation (17) of Li & Ma (1983). The VERITAS observations and analysis results are 
detailed in Table 2. Integral upper limits at 99% confidence level are calculated using the 
Rolke et al. (2005) method, assuming a photon index of Y = 3.0 for the differential power-law 
spectrum dN/dE = N(E/E 0 )~ l . This index was assumed as a moderately softened index as 
compared to the Fermi - LAT index range (1.74—2.9), expected due to absorption by the EBL 
and possible intrinsic turnover. Additionally, this index value is representative of a typical 
TeV blazar. Without detection, the real index in the VHE band remains unknown and 
any integral upper limit derived is dependent on the index assumed. The differential upper 
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limit is quoted at the decorrelation energy, the energy where the calculated flux has minimal 
dependence on the index. Changing the spectral index by ±0.5 changes the differential upper 
limits by less than 10%. 

The results are independently reproduced with two different analysis packages, as de- 
scribed in Cogan (2008) and Daniel (2008). The upper limits range from 1 to 3% of the 
integral Crab Nebula flux above the threshold energy. The energy threshold for each obser- 
vation is defined as the energy at which the differential rate of reconstructed gamma rays 
from the postulated source reaches its maximum and is accurate to within the 20% energy 
resolution of the instrument, a value that is dependent on the observation zenith angle and 
sky brightness. 


4.2. Fermi LAT 

The Fermi - LAT is a pair-conversion telescope sensitive from 20 MeV to > 300 GeV, 
which operates in survey mode. Further details about the characteristics and performance of 
the LAT Instrument on the Fermi Gamma-ray Space Telescope can be found in Atwood et ah 
(2009). Presented here is the analysis of the Fermi LAT data for the six candidates described 
in Section 3. Although the targets were selected based on only 11 months of data, more data 
was available after the completion of VERITAS observations and this larger data set is used 
for the modeling. More specifically, the LAT data from the time period of 2008 August 4 
to 2011 January 4 (MJD 54682.7-55565.0) were used for the modeling analysis. Except for 
variable sources, the analysis procedure was identical for each of the sources and proceeded 
as follows. 

For each candidate, events were extracted from a region of interest (ROI) of 10° radius 
centered on the target coordinates. Events from the “diffuse class” with zenith angle < 100° 
and energy between 300 MeV and 300 GeV were selected. Only data taken during periods 
when the rocking angle of the satellite was < 52° were used to reduce contamination from 
the Earth limb gamma rays, which are produced by cosmic rays interacting with the upper 
atmosphere. The significance and spectral parameters were calculated using an unbinned 
maximum-likelihood method implemented in the LAT Science Tool gtlike 1 (Cash 1979; 
Mattox et al. 1996). A background model was constructed including nearby gamma-ray 
sources and diffuse emission. All sources within 12° of the central source in the second 
Fermi - LAT catalog (2FGL, Nolan et al. (2012)) were included in the model. The spectra 
of known pulsars were modeled by a power laws with exponential cutoffs. As in the 2FGL 


ScienceTools v9r20p0 with the post-launch instrument response function (IRF) P6_V11_DIFFUSE. 
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catalog, a log-parabola function was used for sources with significant spectral curvature. 
Otherwise, spectra were described as a power law. The spectral parameters of the sources 
in the ROI were left free during the fitting procedure. Sources outside the ROI, but within 
the 12° range had their spectral parameters fixed to the 2FGL catalog values. The Galactic 
diffuse emission and an isotropic component, which is the sum of the extragalactic diffuse 
gamma-ray emission and the residual charged particle background, were modeled using the 
recommended hies 5 . 

The LAT data were first analyzed to calculate the time-averaged gamma-ray flux and 
spectral parameters of each candidate. A second analysis was then performed to study the 
impact of the Sun, a bright gamma-ray source, on the flux of candidates located near the 
plane of the ecliptic (RGB J0316+090 and RGB J0909+231). Removing time intervals when 
the Sun was in the ROI of each candidate had a negligible effect on the analysis results. 

Spectral points and a light curve were calculated for each candidate, and a temporal 
analysis was performed to search for flux variability. The timescale of this analysis was 
adjusted based on the specific candidate flux levels. The flux in each energy or time bin was 
determined with the spectral indices of all sources fixed to the best-fit values over the full 
energy and time interval. For an energy or time bin with a test statistic (TS; see Mattox et al. 
(1996)) less than 9 or fewer than 3 predicted photons ( N pred ), a 95% confidence level upper 
limit was calculated. 

The light curves were analyzed to search for flux variability with a likelihood method 
assuming a constant flux for the null hypothesis, following the same procedure as used 
in the 2FGL catalog. Only two sources showed significant evidence of flux variability: 
RGB J0316+090 and RGB J1058+564. For these two sources, a refined analysis was done, 
selecting time periods contemporaneous with the VERITAS observation windows. The dura- 
tion of this contemporaneous period was chosen such that a significant detection (> 5cr) could 
be attained, resulting in slightly extended windows of MJD 55055-55145 for RGB J0316+090 
and MJD 55160-55185 for RGB J1058+564 with respect to the VERITAS observation win- 
dow. For these candidates, a butterfly corresponding to the lcr confidence interval was used 
to represent the spectral information (Figure 3). 

The derived spectral indices for the differential power laws obtained for the candidates 
are relatively hard, ranging from 1.74 to 2.09, as compared to the 2FGL average spectral 
index of 2.21 ± 0.01. The integral fluxes above 300 MeV range from 1.55 to 12.4xl0 _9 ph 


5 The files used were gll_iem_v02_P6_Vll_DIFFUSE . f it for the Galactic diffuse and 

isotropic_iem_v02_P6_Vll_DIFFUSE.txt for the isotropic diffuse component available at 

http : //f ermi . gsf c . nasa . gov/ ssc/data/p6vl 1/ access/lat/BackgroundModels . html 
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cm 2 s 1 , indicating that these six blazars are bright in the high-energy band. The detailed 
Fermi-LAT results for each of these hard-spectrum, bright BL Lacs are summarized in Table 
3. 


4.3. Swift XRT 

The X-ray Telescope (XRT) onboard the Swift satellite (Gehrels et al. 2004) is a focusing 
X-ray telescope sensitive to photons with energy between 0.2 and 10 keV. The data used 
for the broadband SED modeling were analyzed as described in Burrows et al. (2005) with 
the HEASoft package Version 6.9 and XSPEC 6 Version 12.6.0. All data were taken in photon 
counting mode and pile up effects are accounted for when count rates exceeded 0.5 counts 
per second through the use of an annular source region, with a 1-2 pixel inner radius and a 20 
pixel outer radius. Each observation is binned and fit with an absorbed power law between 
0.3 and 10 keV, with the neutral hydrogen density taken from the Leiden/ Argentine/Bonn 
survey of Galactic HI (Kalberla et al. 2005). 

X-ray variability is commonplace for both VHE detected and non-detected blazars. If 
more than one exposure exists for an object and no variability is detected, the de-absorbed 
power-law fit of the combined data set is used to constrain the SED modeling. If variability 
is observed between multiple exposures, results from these separate exposures are shown 
independently on the SED plot, as is the case for RGB J0136+391 and RGB J1058+564. 
These blazars show flux variability factors of ~2 and 3 between exposures, respectively. 
The multiple exposures taken on RGB J0909+231 do not provide sufficient statistics for 
application of an absorbed power-law model and are therefore summed before fitting. Using 
Cash analysis (Cash 1979) did not improve the fitting of the single low-statistics spectra. 
The summed exposure fit result is shown on the SED. 

Each of the absorbed power-law fits applied to the XRT data resulted in photon indices 
greater than 2, with 2-10 keV integral flux levels between 0.3 and 21xl0 _12 ergs cm _2 s _1 . 
The index values suggest that the synchrotron component peaks below keV energies, char- 
acteristic of ISP blazars. The analysis results for each observation are summarized in Table 
4. 


6 http : / /heasarc .nasa.gov/docs/software/llieasoft/xanadu/xspec/XspecManual . pdf 
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4.4. Swift UVOT 

The Swift-X . RT observations were supplemented with simultaneous UVOT exposures 
taken in the V, B, U, UVW1, UVM2, and UVW2 bands (Poole et al. 2008). The UVOT 
photometry is performed using the HEASoft program uvotsource. The circular source region 
has a 5" radius and the background region consists of several 15" radii circles of nearby 
empty sky. The results are reddening corrected using E(B-V) coefficients (Schlegel et al. 
1998). The Galactic extinction coefficients are applied according to Fitzpatrick (1999). The 
uncertainty in the reddening E(B-V) is the largest source of error, especially in the UV bands 
for blazars that have a large value of E(B-V). If more than one exposure exists in a specific 
band for an object, the data from the observation closest to the VERITAS exposure are 
used, although no significant variability is seen across any band for any blazar. A summary 
of the UVOT analysis results is presented in Table 5. 


5. Broadband SSC Modeling 

Leptonic models for blazar jet emission attribute the higher energy peak in the SED 
to the inverse-Compton scattering of lower energy photons off a population of non-thermal, 
relativistic electrons. These same electrons are responsible for the lower-energy synchrotron 
emission that makes up the first peak. The target photon field involved in the Compton 
upscattering can either be the synchrotron photons themselves, as is the case in SSC models, 
or a photon field external to the jet in the case for EC models. 

The previously described multiwavelength data are matched with archival radio data 
collected from NED. These data are used to test a steady-state leptonic jet model for the 
broad-band continuum emission from the blazars. Although it has been found that ISP 
BL Lacs are sometimes better represented by external- Compton models, e.g. Acciari et al. 
(2009b), taking into account the lack of redshift information and lack of constraints from the 
broadband data we prefer not to apply an EC model, which would introduce additional free 
parameters, compared to the SSC model applied here. The model-predicted flux reflects the 
absorption of VHE gamma rays by the EBL according to the redshift information summarized 
in Table 1 using the model from Gilmore et al. (2009). The level of TeV absorption resulting 
from this model is consistent with the absorption derived from the Finke et al. (2010) and 
Franceschini et al. (2008) models. 

The SSC model applied to the multiwavelength data is the equilibrium version of the 
Bottcher & Chiang (2002) model, as described in Acciari et al. (2009b). In this model, the 
emission originates from a spherical blob of relativistic electrons with radius R. This blob is 
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moving down the jet with a Lorentz factor F, which corresponds to a speed of /3rc. The jet 
axis is aligned toward the line of sight with an angle 9 0 b s , which results in Doppler boosting 
with a Doppler factor D = (T[l — fir cos 0 o bs]) _1 - In order to minimize the number of free 
parameters, we assume that 9 obs = 1/T, often referred to as the critical angle, for which 
T = D. 

Within the model, nonthermal electrons are injected and accelerated into a power-law 
distribution Q( 7) = <5o7 _<? between the low- and high-energy cut-offs, 7 m i n and 7 max - The 
radiation mechanisms considered lead to an equilibrium between particle injection, radiative 
cooling and particle escape. This particle escape is characterized with an escape efficiency 
factor 77, such that the time scale of escape t e sc = r/Rjc. This results in a particle distribution 
which streams along the jet with a power L e . Synchrotron emission results from the presence 
of a tangled magnetic field B, with a Poynting flux luminosity of L B - 

The two parameters L e and L B allow the calculation of the equipartition parameter 
e Be = L B /L e . This equipartition parameter is used as an estimator of the feasibility of 
the model, where models which result in e Be ~ 1 are preferred. If the particle energy 
density greatly dominates over the magnetic field energy density, namely a particle dominated 
scenario, then the magnetic field cannot serve to collimate the jet. Following this design, 
acceptable parameters should result in at least partition conditions with L B > L e . 

The broadband SED for each blazar can be seen in Figure 3, with the SSC model 
parameters for each representation summarized in Table 6. For each blazar, the archival 
radio data are taken as upper limits as these measurements are believed to contain a large 
amount of radiation produced in the radio lobes in addition to the synchrotron emission from 
the jet. The modeling for each of the six blazars shows synchrotron peak locations ~ 10 15 
Hz, characteristic of borderline ISP/HSP blazars. 

RGB J0136+391: This blazar is modeled using the lower of the two variable X-ray 
states for three different assumed redshifts (z = 0.2, 0.3 and 0.4). The variability timescale of 
three months suggested by the factor-of-two variability observed between XRT exposures is 
not short enough to provide a constraint on the size of the model emission region. Under the 
assumption that the gamma-ray emission remained steady during the Fermi and VERITAS 
observations, the hard LAT spectrum and low VERITAS upper limit derived for this blazar 
suggest a steepening of the gamma-ray spectrum at E > 100 GeV, which could be caused 
by the EBL absorption if a redshift * > 0.4 is assumed. Alternatively, this apparent break 
could originate from uncorrelated variability in the high and very-high energy bands. More 
specifically, a low flux state in the VHE band during VERITAS observations could provide 
a redshift- independent explanation of the apparent spectral softening. Figure 3 shows the 
model predictions, corrected by EBL absorption, for each of the redshift values. Only the 
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model at z = 0.4 is compatible with the VERITAS upper limit, which also results in a 
framework with balanced radiation and particle energy. 

RGB J0316+090: Due to the variability detected in the high-energy band, this blazar 
is modeled with Fermi - LAT data which spans the complete time period sampled by shorter 
VERITAS observations. The UVOT errors resulting from the E(B-V) reddening correction 
for this blazar are so large that all exposures except the V band are unconstraining to the 
SED modeling and therefore not shown. Application of the SSC model for an assumed 
redshift of z — 0.2 results in a particle dominated framework consistent with the VERITAS 
upper limit. 

RGB J0909+231: The SSC model parameters used to describe the broadband emis- 
sion of this blazar are determined for a redshift of z = 0.5, based on the lower limit derived 
from Mg II absorption lines found in the SDSS data. The Swift X-ray data from three 
exposures have been summed in order to provide sufficient statistics for absorbed power-law 
fitting. The model results in a particle-dominated scenario with a slight discrepancy between 
the model and the Fermi upper limit in the 1 — 3 GeV energy bin. Moving the Compton 
peak accomodates this upper limit only results in an even less favored, particle-dominated 
emission state. 

RGB J1058+564: Due to the variability observed in the LAT band, the SED modeling 
for this blazar is done with LAT data from the time period coincident with the VERITAS 
observation window and for the lower of two X-ray states observed. The variability timescale 
suggested by the variability factor of three observed between XRT exposures is not short 
enough to provide any constraint on the size of the emitting region. Both X-ray states are 
shown on the SED for reference. The SSC model shows agreement with the broadband data, 
although results in a disfavored particle dominated scenario. 

RGB J1243+364: The modeling for this blazar is completed for a redshift of z = 0.5, 
based on the new lower limit found from Mg II absorption lines in the public SDSS data. 
The model agrees with the broadband data and permits parameters at equipartition. 

RX J 1436.9+5639: The modeling for the broadband data of this blazar allows marginal 
agreement with the Fermi data when completed for the redshift of the spatially coincident 
galactic supercluster, more specifically z = 0.15. The resulting model parameters are far 
below equipartition, suggesting a particle dominated jet. 
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6. Discussion 

Six promising TeV blazar candidates were selected from the 1FGL catalog for observa- 
tions with VERITAS. These observations resulted in no VHE detections. Initial selection 
of these candidates from 1FGL power-law fit extrapolation suggested integral flux levels be- 
tween 1 and 12% of the Crab Nebula flux above 150 GeV after accounting for absorption 
by the EBL. With additional Fermi - LAT data and more information about the blazar red- 
shifts, the expected fluxes were updated to levels between 0.3% and 3.5% of the Crab Nebula 
flux above 150 GeV. The VERITAS exposure times were allocated based on the initial ex- 
trapolated values, resulting in only one upper limit below the updated VHE extrapolation, 
namely the upper limit for RGB J0136+391. The non-detection of this blazar suggests spec- 
tral steepening of the high-energy spectrum that can either be explained by intrinsic spectral 
curvature, redshift-dependent EBL absorption or uncorrelated variability. 

Multiband observations are presented and allow for the construction and SSC modeling 
of the radio through TeV broadband SEDs. The model applied to these blazars is suffi- 
cient to represent the broadband data for each of the six, with model parameters roughly 
comparable to those found for other VHE-detected ISPs using the same model (Abdo et al. 
2011; Acciari et al. 2009b). The magnetic fields obtained in the modeling are generally low, 
resulting in disfavored, particle-dominated jets. This condition could be relaxed by including 
an external photon held for inverse-Compton scattering, allowing solutions closer to equipar- 
tition, as was done in Abdo et al. (2011) and Acciari et al. (2009b). However, this scenario 
has not been explicitly tested. The quality of the data sampling of the inverse-Compton 
peak and the fact that the redshifts are not well determined for the majority of the objects 
do not provide sufficient constraints for a model with the extra degrees of freedom associ- 
ated with adding an external source of seed photons. Similar parameters were also found 
for the borderline ISP/HSP TeV detected blazar PKS 1424+240, with the exception of the 
spectral index for the injected electron distribution, which was found to be very soft ( q = 
5.1, Acciari et al. (2010)), while the typical values for this model fall between q =2.3 and 
2.7. 


Each of these blazars shows a synchrotron peak frequency characteristic of the ISP/HSP 
divide, namely u synch 10 15 . Additionally, each of these BL Lacs exhibits a similar flux 
level within the high-energy gamma-ray band, showing comparable L sync h and L ic- This 
commonality of sub-type and peak luminosities is likely a bias introduced to the selection 
process through the sensitivity of the Fermi LAT instrument being greater in the 300 MeV 
-100 GeV band as compared to the >100 GeV band. Nearby blazars releasing a large 
fraction of power above 100 GeV are known to be good candidates for detection in the VHE 
regime. This type of power emission is characteristic of HSP BL Lacs, the most commonly 


18 


VHE-detected type of blazar, but the least frequently detected BL Lac in the 1FGL and 
2FGL catalogs. 

Although the selection of these six blazars from high-energy Fermi - LAT data did not 
lead to any new TeV blazar discoveries, the exercise has been very illuminating. We 
are reminded that when selecting candidates for observation with TeV instruments, ac- 
curate redshift information is crucial. Without this information it is difficult, to deci- 
pher if the lack of TeV emission is due to the intrinsic emission mechanism or the ab- 
sorption of gamma rays by the EBL. The redshift lower limits that were found in the 
SDSS for RGB J0909+231 and RGB J1243+364 indicate distances where future TeV de- 
tection from a deeper exposure is unlikely. Continued TeV observations of RGB J0136+361, 
RGB J0316+090 and RX J1436. 9+5639, the three sources that remain without solid red- 
shift information, could provide insight into the distance to these sources, while further 
observations of RBG J1058+564 can be directly applied to further investigate the emission 
mechanism at work within the blazar jet. 
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Fig. 1. — Optical spectra for the six BL Lacs selected from the 1FGL catalog and observed 
with VERITAS. The black shows the object spectrum, while the red shows the instrumental 
noise. Only one BL Lac had a confirmed redshift upon selection (RGB J1058+564, z=0.143), 
confirmed with the SDSS spectrum shown in the middle-right panel. Redshift lower limits for 
RGB J0909+231 [z > 0.4305; middle left) and RGB J1243+364 (z > 0.485; lower left) are 
found from Mg II absorption lines in the SDSS spectra. A featureless SDSS spectrum is found 
for RX J1436. 9+5639. A redshift measurement attempt for RGB J0136+391 using the Keck 
LRIS instrument shows a featureless power- law spectrum (upper left). The normalized ESI 
spectrum of RGB J0316+090 exhibits three unidentified absorption features (upper right). 
The spectral analysis of the Keck LRIS and ESI spectra are detailed in Kaplan et al. (2012). 
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Fig. 2. — Fermi - LAT light curves, with units of 10 -9 ph cm -2 s -1 , are shown for the six can- 
didate VHE-emitting BL Lacs for 29 months of LAT data (MJD 54682-55565; 2008 August 
4-2011 January 4). The beginning and end of the VERITAS observations are denoted by ver- 
tical grey lines in each panel. The short VERITAS observation periods for RGB J0136+364 
and RGB J0316+090 can be seen, representing only 7 days each. Upper limits at 95% confi- 
dence level are shown for bins resulting in TS of less than 9, denoted by downward pointing 
black triangles. Only RGB J0316+090 and RGB J1058+564 show any significant variability. 
For these two sources, LAT data only for the time periods within the window of VERITAS 
observations are used to constrain the modeling. 
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Fig. 3. — Non-contemporaneous broadband SED data for each BL Lac shown with corre- 
sponding SSC modeling using the model of Bottcher & Chiang (2002). The modeling results, 
corrected for the EBL absorption, are shown with grey lines. The archival radio data points 
are taken from NED and used only as upper limits, accounting for the fact that much of the 
radio emission may result from diffuse synchrotron emission in the radio lobes of the jet. See 
the text for a more detailed description of the model parameters, with values listed in Table 
5. 
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Table 3. Summary of Fermz-LAT observations and analysis results. A variability study was completed using 29 
months of data for all sources. For steady sources, the spectral analysis is completed for 29 months of data. The two 
variable sources RGB J0316+090 and RGB J1058+564 show spectral analysis results for LAT data coincident with a 
slightly expanded window as compared to the VERITAS observation window so as to allow a 5cr detection. 
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Table 5. Summary of Swift - UVOT observations and analysis results. 


Target 

Observation 

ID 

Date 

[MJD] 

Band 

Frequency 

[Hz] 

vF u 
[Jy Hz] 

Used in 
SED? 

RGB J0136+391 

00039107001 

55134 

UVWl 

1.14 x 10 15 

(1.16 + 0.15) x 10 12 

y 


00039107001 

55134 

UVM2 

1.34 x 10 15 

(1.44 + 0.30) x 10 12 

y 

RGB J 0316+090 

00038370001 

54899 

V 

5.55 x 10 14 

(2.24 + 1.04) x 10 12 

y 


00038370001 

54899 

U 

8.57 x 10 14 

(2.92 + 3.87) x 10 12 

n* 


00038370001 

54899 

UVWl 

1.14 x 10 15 

(2.44 + 7.64) x 10 12 

n* 


00038370001 

54899 

UVW2 

1.48 x 10 15 

(2.8 + 18.5) x 10 12 

n* 

RGB J0909+231 

00040540001 

55338 

V 

5.55 x 10 14 

(2.35 + 0.11) x 10 11 

y 


00040540001 

55338 

B 

6.93 x 10 14 

(1.77 + 0.23) x 10 11 

y 


00040540001 

55338 

U 

8.57 x 10 14 

(2.26 + 0.21) x 10 11 

y 


00040540001 

55338 

UVWl 

1.14 x 10 15 

(1.78 + 0.18) x 10 11 

y 


00040540001 

55338 

UVM2 

1.34 x 10 15 

(1.63 + 0.22) x 10 11 

y 


00040540001 

55338 

UVW2 

1.48 x 10 15 

(1.80 + 0.19) x 10 11 

y 


00040540002 

55338 

UVM2 

1.34 x 10 15 

(1.90 + 0.22) x 10 11 

n 


00040540003 

55338 

UVW2 

1.48 x 10 15 

(1.74 + 0.20) x 10 11 

n 

RGB J 1058+564 

00038215001 

54852 

V 

5.55 x 10 14 

(1.42 + 0.04) x 10 12 

y 


00038215001 

54852 

U 

8.57 x 10 14 

(1.45 + 0.04) x 10 12 

y 


00038215001 

54852 

UVW2 

1.48 x 10 15 

(1.49 + 0.05) x 10 12 

y 


00038453001 

54888 

B 

6.93 x 10 14 

(1.64 + 0.07) x 10 12 

n 


00038453001 

54888 

U 

8.57 x 10 14 

(1.70 + 0.06) x 10 12 

n 


00038453001 

54888 

UVWl 

1.14 x 10 15 

(1.66 + 0.06) x 10 12 

n 

RGB J1243+364 

00038445001 

54875 

V 

5.55 x 10 14 

(1.28 + 0.05) x 10 12 

y 


00038445001 

54875 

B 

6.93 x 10 14 

(1.30 + 0.04) x 10 12 

y 


00038445001 

54875 

u 

8.57 x 10 14 

(1.42 + 0.05) x 10 12 

y 


00038445001 

54875 

UVWl 

1.14 x 10 15 

(1.34 + 0.05) x 10 12 

y 


00038445001 

54875 

UVM2 

1.34 x 10 15 

(1.49 + 0.05) x 10 12 

y 


00038445001 

54875 

UVW2 

1.48 x 10 15 

(1.47 + 0.05) x 10 12 

y 

RX J1436. 9+5639 

00038435001 

54918 

V 

5.55 x 10 14 

(9.03 + 0.09) x 10 11 

y 


00038435001 

54918 

B 

6.93 x 10 14 

(1.14 + 0.07) x 10 12 

y 


00038435001 

54918 

U 

8.57 x 10 14 

(9.71 + 0.05) x 10 11 

y 


00038435001 

54918 

UVWl 

1.14 x 10 1S 

(1.06 + 0.06) x 10 12 

y 


00038435001 

54918 

UVM2 

1.34 x 10 15 

(1.21 + 0.07) x 10 12 

y 


00038435001 

54918 

UVW2 

1.48 x 10 15 

(1.12 + 0.06) x 10 12 

y 


00038289001 

54918 

V 

5.55 x 10 14 

(8.58 + 0.05) x 10 11 

n 


00038289001 

54918 

U 

8.57 x 10 14 

(1.06 + 0.04) x 10 12 

n 


00038289001 

54918 

UVW2 

1.48 x 10 15 

(1.26 + 0.06) x 10 12 

n 


^Galactic reddening dominates the uncertainty with E(B-V) = 0.356 for RGB J0316+090; U, UVW1, 
UVW2 band flux measurements do not constrain the SED. 


Table 6 . Summary of the broadband SED SSC modeling parameters. See text for parameter descriptions. 
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